Chaos and synchronized chaos in an earthquake model 
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We show that chaos is present in the symmetric two-block 
Burridge-Knopoff model for earthquakes. This is in contrast 
with previous numerical studies, but in agreement with exper- 
imental results. In this system, we have found a rich dynam- 
ical behavior with an unusual route to chaos. In the three- 
block system, we see the appearance of synchronized chaos, 
showing that this concept can have potential applications in 
the field of seismology. 



In recent years, the phenomenon of chaotic synchro- 
nization has been a subject of intensive studies. By def- 
inition, chaotic systems present strong sensitivity to the 
initial conditions, and in principle it seems impossible to 
synchronize them. However, Fujisaka et al. [Q and Pecora 
et al. Q showed that systems with chaotic behavior can 
be synchronized, if appropriate connections among them 
are made. This phenomenon has been called "chaotic 
synchronization" , and has been investigated mainly in 
applications for secure communications ||. 

Another area of active research nowadays is related 
to systems that present avalanche-like dynamics. This 
was triggered by a paper by Bak, Tang and Wiesenfeld 
||. They showed that certain dissipative systems, with 
many degrees of freedom, naturally evolve to a critical 
state characterized by power-law distributions in space 
and time. They denoted this phenomenon Self-Organized 
Criticality (SOC). 

One of the systems that have been studied in con- 
nection with SOC is the Burridge-Knopoff (B-K) model 
for earthquakes f|. This model consists of blocks con- 
nected by springs. The whole systems is pulled with con- 
stant velocity on a surface with friction. It has been 
shown experimentally |5j and numerically that this 
model presents a region of power-law distribution simi- 
lar to what is observed in real earthquakes, namely, the 
Gutenberg- Richter law . Since the power-law does not 
span the entire system, one could say that this system 
does not present what has been denned as SOC. How- 
ever, a variation of it, called the "train model", does 
present SOC §,§. 

After the work by Carlson and Langer , several stud- 
ies on the B-K model were performed. With respect to 
the chaotic properties of the model, we are aware of the 
numerical studies by Nussbaum and Ruina |n|, Huang 
and Turcotte pdj| , Nakanishi p^ ], Crisanti et al. and 
Lacorata et al. |l4f| . In jl2| and Jl3| the systems consid- 
ered had more than two blocks and they were evolved 



by cellular automaton rules. Nussbaum et al. studied a 
symmetric two-block B-K model, and verified that, if a 
friction force of the Coulomb type (that is the dynamic 
friction coefficient being constant) the system presents 
only period behavior. Huang and Turcotte, and Lacorata 
and Paladin found chaotic behavior in the two-block B-K 
model only with the presence of an asymmetry in the sys- 
tem, even considering a velocity weakening friction force 
In particular, they considered the friction force 
in one block being different from the friction force in the 
other block. On the other hand, by modeling the two- 
block B-K model by electronic circuits, Field, Venturi 
and Nori showed experimentally that a completely 
symmetric system does present chaotic behavior in a wide 
range of the parameter space. Therefore, their results are 
in contradiction to what was reported in jn]|y|]. One of 
the purposes of this letter is to resolve this contradiction. 
We show that the two-block B-K system in a symmetric 
configuration is chaotic. As in the experimental study, 
chaos is seen in wide range of parameter values. We 
stress that a one-block system in the BK model (with a 
linear drive) cannot present chaos, since its dimension- 
ality is smaller than the minimum dimension necessary 
for a system to present chaotic behavior, which is three. 
(It is obvious that here we are considering the absence of 
elements with delay in the system, since in this way one 
could increase its dimensionality up to infinity). 

In this letter we also show that the phenomenon of 
synchronized chaos appears in the three-block system of 
the Burridge-Knopoff model for earthquakes. Most im- 
portantly, it comes naturally from the geometry of the 
system, without any need for special connections, as it is 
generally the case in the studies of synchronized chaos. 
From our results, we speculate that synchronized chaos 
may have applications in the field of seismology. That is, 
earthquakes faults, which are generally coupled through 
the elastic media in the Earth crust could in principle 
synchronize, even when they have an irregular (chaotic) 
dynamics. As a consequence of synchronization, the di- 
mensionality of the system decreases what simplifies the 
analysis of the system, as discussed below. We quantify 
the degree of synchronization by studying the Liapunov 
exponent associated with the synchronization manifold, 
that is, the transverse (or conditional) Liapunov expo- 
nent PlJiql and compare it with the Liapunov exponents 
of the three-block system. 

We start by reviewing the B-K model. It consists of 
a chain of blocks of mass m, connected by coil springs 
of strength k c to its nearest neighbors. They are situ- 
ated on a rough surface. Between the blocks and the 
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surface there is a frictional force F. Here we consider 
that F is a function of the block velocity. Each block 
is also attached by a leaf spring of strength k p to a 
surface that moves with constant velocity v. A figure 
of the system can be found in |(|. Following Carlson 
and Langer Q , we denote by Xj the position of block j 
with respect to its equilibrium position, write the friction 
force as F(Xj/v c ) — F 3>(Xj/v c ) where <J>(0) = 1 and v c 
is a characteristic velocity, and introduce the variables 
r = ujpt, ujp = kp/m, Uj = k p Xj/F . Then, the equation 
of motion for the two-block system can be written in the 
following dimensionless form 
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with v = v/V , v c = v c /V , V = F J y/k p m and 
k = k c /k p . Dots denote differentiation with respect to 
r. Eq. [l] is valid only when block j is moving. If it is 
not, its equation is simply Uj = v. We use the velocity 
weakening friction force introduced in j6j, given by 



1 



l + U/v c 



(2) 



which is a simple nonlinear function. In the simulations 
displayed here, we did not allow backward motion, that is 
the static friction force can take any value to prevent it. 
However, several tests showed that if backward motion is 
allowed, the results remain essentially the same. 

In contradiction with the reports in but in 

agreement with experimental results |l5| , we find that 
a velocity weakening friction force leads to a rich dynam- 
ical behavior in the two-block B-K system, even when it 
is identical for the two blocks. We find the presence ofpe- 
riodic, quasiperiodic and chaotic behavior. In Fig. |lj(a) 
we show an example of chaotic motion, by plotting U\ 
versus U t - Uf, with l/v c = = 1/ff = 1. We de- 
note Uj = vt — v c /(i> c + v) the unstable equilibrium point 
around which the orbits of block j circle in phase space, 
which is found by taking Uj = and Uj = v in Eq. ^ 
(the stability of such a solution for any number of blocks 
was analyzed in ||). Unless explicitly stated otherwise, 
we take in the numerical studies shown here k = 1 and 
v = 0.1. However, similar behavior was found for other 
parameter values, as well. We do not display U 2 versus 
U 2 — U 2 , since we have found that for this system the at- 
tractors of the two blocks are the same. In other words, 
the plot of U 2 versus U 2 — U 2 is identical to the one of 
Ui versus U\ — f7f , showing that the two blocks will visit 
the same region of the phase-space, but not necessarily at 
the same time. The initial conditions for the simulations 
shown here arc the blocks initially at rest and with small 
random displacements from their equilibrium position. 
In all the simulations, a transient time was discarded. 

The bifurcation diagram for the two-block system is 
shown in Fig. p](b). There, we display U\ versus l/v c , 



in the Poincare section satisfying U\ — Uf = 0. In or- 
der to quantitatively characterize the dynamics, we have 
calculated the two largest Liapunov exponents of the sys- 
tem. If the largest Liapunov exponent (LLE) of the sys- 
tem is greater than zero, then, by definition the system 
is chaotic. Quasiperiodic motion occurs in this system 
when the LLE and the second largest Liapunov expo- 
nent (SLLE) are zero, and the motion is periodic when 
the LLE is zero and the SLLE is negative. We used 
the method introduced in jl9) to calculate the LLE and 
SLLE, and our results are displayed in Fig. |l](c). The 
solid line refers to the LLE and the dashed line is the 
SLLE. We investigated in detail the first entrance into 
chaos, which occurs at \ jv c w 0.112 for these parameter 
values, and noticed that it is unusual. We found that the 
route to chaos is from period one, to period two and then 
directly into chaos. This unusual route is probably due 
to the fact that here we have a system governed by non 
differentiable flows. Most of what is known in dynami- 
cal system theory deals with systems that are infinitely 
differentiable, which is not the case here. More details 
about return maps and other quantities for this system 
will be published elsewhere |T^ ]. 

One may ask, why Huang and Turcotte jll]] and Laco- 
rata and Paladin |Q did not find chaotic behavior in the 
symmetric B-K model? The reason is that they made an 
inconsistent assumption, that is, that the driving block 
does not move during the slipping events. This assump- 
tion is equivalent to take the pulling velocity going to 
zero, since they drop out the term vt in Eq. [ll However, 
from the equations of motion of the system, it turns out 
that the blocks will not move, if they are initially at rest, 
if one considers this, since the displacement of the blocks 
is proportional to the pulling velocity. This is shown in 
detail in ||. On the other hand, the blocks will not stop 
if they are initially in motion. That is, stick-slip dynam- 
ics does not occur when the system is evolved by the 
equation used by those authors. One way to avoid this is 
to take a discontinuous friction force as in |l8[ , and this 
was not considered in the study of the symmetric system 
by Huang and Turcotte, and Lacorata and Paladin. In 
other words, their equations of motion does not describe 
the dynamics of the Burridge-Knopoff model. 

Now, we concentrate on the phenomenon of chaotic 
synchronization in the B-K model. We have not found 
this phenomenon in the two-block system with the fric- 
tion force we consider here. If the two blocks were syn- 
chronized, they would behave as a single block, and this 
would be equivalent to have a single block system with 
chaotic behavior, which we know it is impossible, as dis- 
cussed above. We see however, that a three block system 
does present chaotic synchronization. 

The equation of motion for the three-block system in 
dimensionless form is 



Ui = k(U 2 - Ui) -Ui + 
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U 2 = k{Ux - 2U 2 + U 3 ) -U2 + VT- $>{U 2 /v c 2 ), 

Us = k(U 2 - U 3 ) -U 3 + vt- $(U a /v®. (3) 

We see in the equations that govern the motion of blocks 
1 and 3 have the same functional form. They are also 
linked to a common subsystem, that is, block 2. This con- 
figuration is not of the "master-slave" type, since there 
is feedback between blocks 1 and 2 and between blocks 2 
and 3. However, the equations have the necessary ingre- 
dients for chaotic synchronization between blocks 1 and 
3 to occur, which is the same functional form. 

We have found chaotic synchronization when the pa- 
rameters for all the blocks are the same only in a very 
small region of the parameter space. It is not difficult to 
understand why. With the leaf springs having the same 
value, blocks 1 and 3 are more loose than block 2 (since 
the later is attached to two coil springs instead of one). 
Therefore, in general, blocks 1 and 3 attain larger veloci- 
ties than block 2, and with the friction force that we use, 
they are more unstable. Chaotic synchronization hap- 
pens when the subsystems to be synchronized are be more 
stable than the subsystem to which they are connected 
. Also, it is a property of chaotic synchronization that 
blocks 1 and 3 need to have the same parameter values, 
if perfect synchronization in the absence of control is the 
goal. 

However, we find that when there is an asymmetry 
in the system, chaotic synchronization occurs in a large 
range of the parameter space. We can either introduce 
asymmetries in the springs, masses, or friction forces. We 
choose the later to demonstrate our results, but similar 
results were found in the other two cases. Therefore, we 
make the friction force in block 2 smaller than the friction 
force in blocks 1 and 3 (which means that they are more 
rough than block 2). 

In Fig. ||(a) we show an example of a chaotic orbit 
for block 1, with the parameter values \/v c = \jv\ = 
^l v 3 = 4/^2 = 0.165. We will see that this attractor is 
from a regime in which blocks 1 and 3 are synchronized. 
In Fig. |(b) the LLE (solid line) and the SLLE (dashed 
line) of the three-block system is shown. In Fig. ||(c) 
we show the Liapunov exponent of the synchronization 
manifold (solid line), that is, what is called the transverse 
(or conditional) Liapunov exponent p]Jlq|. We have cal- 
culated the transverse Liapunov exponent by adapting 
the method introduced by Benettin ]19| . That is, af- 
ter the transient dies out, we evolve the orbit of block 
3 by making it slightly different from the orbit of block 
1. Then, we verify how the difference between the orbits 
of the two blocks evolves after a short time step. The 
perturbation is renormalized in the direction of the max- 
imum growth, and the process is repeated many times. 
The transverse Liapunov exponent is given by the av- 
erage logarithm (in this paper we use base 2) of the 
growth (or shrinkage) of the perturbation along the or- 
bit. Fig. ||(c) also shows (dashed line) the Euclidean dis- 



tance D in phase-space between blocks 1 and 3, that is, 

D = yVi - U 3 ) 2 + (t/i - t/ 3 ) 2 , as a function of l/u°. 
This distance is an average over a time At = 2000. We 
see that the transverse Liapunov exponent correctly de- 
termines the region in which blocks 1 and 3 are synchro- 
nized. In this situation, the transverse Liapunov expo- 
nent is negative and the distance between the two blocks 
is zero. Comparing Fig. |2|(b) and Fig. ||(c) we can identify 
the regions of chaotic synchronization, where we have one 
or more Liapunov exponents of the three-block system 
greater than zero and the transverse Liapunov exponent 
is less than zero. This is, for example, the case of the 
attractor shown in Fig. ^](a). There, only one Liapunov 
exponent of the system is larger than zero. The com- 
parison between Fig. ||(b) and Fig. ||(c) also shows that 
the transition from chaos to hyperchaos |^(| does not de- 
termine here the transition from chaotic synchronization 
to non-synchronization, as it was the case of the system 
studied in (2^] (the hyperchaos regime is defined as the 
one in which there is more than one positive Liapunov 
exponent). For example, at l/v c = 0.38 there is only 
one positive Liapunov exponent, and no synchronization 
is seen between blocks 1 and 3. There is also the case 
in which the LLE and the SLLE are non-positive, and 
blocks 1 and 3 do not synchronize, as in \/v c = 0.4. 

Now we discuss what is the possible relevance of our 
finds to the analysis of real earthquakes. It is known 
that earthquakes are not strictly periodic phenomena. If 
their irregular behavior is caused by a chaotic, determin- 
istic process, as this simple spring-block model suggests, 
then in principle, prediction about them could be made 
in a short time scale. In this way, the damage they cause 
could be minimized. The problem with modeling chaotic 
systems is that, if the dimension of their attractor is very 
large, very little can be done with respect to prediction 
p2]| . Such systems are treated basically as stochastic 
systems. However, if synchronization happens among 
elements of a large system, then the dimensionality of 
the attractor decreases. For example, in the three-block 
system studied here, the dimension of the attractor can 
decrease by up to two. If the dimension of the attrac- 
tor decreases, the analysis of the system becomes easier. 
This fact emphasizes the need for a deeper analysis of the 
relationship between synchronized chaos and prediction. 
We hope that this work will motivate more studies on 
chaotic synchronization for applications to seismology. 
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FIG. 2. Orbits in phase-space for Ui versus U\ — vt when 
l/v c = l/v{ = \/vg = 4/i;| = 0.165. (b) The largest Lia- 
punov exponent Ai (solid) and the second largest Liapunov 
exponent A2 (dashed) as a function of l/v c . (c) The trans- 
verse Liapunov exponent At (solid line) and the Euclidean 
distance D in phase-space between blocks 1 and 3 (dashed 
line) a function of \jv c . Here we have v = 0.1, fc = 1 in 
a three-block system. The Liapunov exponents were calcu- 
lated with perturbations of 10 -4 , time step of 0.05, and with 
400000 iterations. The dotted line at D — is just a guide to 
the eye. 



FIG. 1. Orbit in phase-space for block 1 when 

ljv\ — l/fl = 1. (b) Bifurcation diagram of U\ in the 
Poincare section in which Ui — U{ = 0. (c) The largest Li- 
apunov exponent Ai (solid line) and the second largest Lia- 
punov exponent A2 (dashed line) of the system. Here we have 
v — 0.1 and k = 1 in a two-block system. The Liapunov 
exponents were calculated with perturbations of 10~ 4 , time 
step of 0.05, and with 400000 iterations. The dotted line at 
A = is just a guide to the eye. 
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